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The dynamics of the annihilation of a vortex- antivortex pair is investigated. The pair is activated 
magnetically during the run of a simulated hysteresis loop on a square mesoscopic superconducting 
cylinder with an antidot inserted at its center. We study the nucleation of vortices and ant ivort ices 
by first increasing the magnetic field, applied parallel to the axis of the sample, from zero until the 
first vortex is created. A further increase of the field pulls the vortex in, until it reaches the antidot. 
As the polarity of the field is reversed, an antivortex enters the scene, travels toward the center of the 
sample and eventually the pair is annihilated. Depending on the sample size, its temperature, and 
Ginzburg-Landau parameter, the vortex-antivortex encounter takes place at the antidot or at the 
superconducting sea around it. The position and velocity of the vortex and antivortex singularities 
were evaluated as a function of time. The current density, magnetization and order parameter 
topology were also calculated. 



Achieving a deep understanding of the nucleation and 
propagation of vortices in real superconductors is a truly 
complex task, since these entities interact with almost 
everything: first, with the surface of the specimen, to 
surpass it; upon entrance, with other vortices that might 
have already penetrated, and also with defects, which 
might attract them and even act as pinning centers. Ad- 
ditional difficulties to emulate the problem arise from the 
fact that vortices generate heat while propagating, what 
can be harmful to the robustness of the superconduct- 
ing properties, if not catastrophic, as is the case of vor- 
tex avalanches observed in some superconducting films 
01 0- S Hi & @- It is quite common, however, that 
the existence of pinning potentials represent a beneficial 
feature, since vortices can thus be prevented from under- 
going dissipative motion. An interesting approach to the 
problem, which enables one to address most specificities 
without excessive complexity, is to work in the small uni- 
verse of mesoscopic samples. In such an ambient, one can 
accommodate the essential ingredients: relatively impor- 
tant surface-to- volume ratio, only a few vortices on scene, 
and a number of defects - the so-called antidots - usu- 
ally arranged in a regular pattern. Furthermore, one can 
study the interaction of an individual vortex-antivortex 
(V-AV) pair and, eventually, witness their mutual anni- 
hilation. 

Recently, there has been many studies about V-AV 
configurations in mesoscopic superconductors (see for in- 
stance References [7|, [8|, 0, [HI E|). The authors of these 
references have found that vortices and antivortices may 
coexist in equilibrium in configurations which look like 
a V-AV molecule. A somewhat common approach is to 
assume an a priori configuration and minimize the free 
energy in terms of some relevant parameter for which 



the V-AV molecule is a stable configuration. Here, we 
will focus in a rather different approach concerning more 
with the dynamics of a V-AV encounter. The aim of the 
present work is to elucidate the details involved in the 
process of creation of pairs, following their time-evolution 
and ultimate disappearance. We opted to do this making 
no use of a priori assumptions regarding symmetries, but 
simply varying the applied magnetic field, allowing for 
the spontaneous nucleation of a vortex and, cycling the 
field, of an antivortex, which form a pair of interacting 
entities, whose subsequent time-evolution is monitored. 

In this Letter we present the dynamics of the anni- 
hilation of a V-AV pair activated magnetically during 
a short hysteresis loop. The geometry we consider is a 
square mesoscopic superconducting cylinder with an an- 
tidot placed at its center. Simulations are made in the 
presence of an applied magnetic field parallel to the axis 
of the cylinder. We study the nucleation of vortices and 
antivortices by first increasing the magnetic field from 
zero until the first vortex is created. The field was then 
decreased toward negative values. During this process 
the vortex travels heading the center of the sample, is 
trapped by the antidot and is then annihilated by an an- 
tivortex as the polarity of the applied field is reversed. 
To monitor the whole process, we evaluate the position 
and velocity of the vortex and antivortex singularities 
as a function of time. We also calculate the order pa- 
rameter topology, the current density, and the magne- 
tization. In order to follow the dynamics of the physi- 
cal quantities we use the gauge- invariant time-dependent 
Ginzburg-Landau (TDGL) equations. These equations 
describe the time-evolution of the complex order param- 
eter ip and the vector potential A, which is related to the 
local magnetic field through the expression h = V x A. 
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In the absence of transport currents, these equations read 

^ = -HV-A) 2 ^ + (l-T)#-H 2 ), 
8 A 

{3— = J s -« 2 Vxh, (1) 

where the supercurrent density is 

J a = (l-r)8[^HV-A)^] . (2) 

Here, the distances are measured in units of the coher- 
ence length at zero temperature £(0); the magnetic field 
in units of the upper critical field at zero temperature 
H c2 (0); the temperature in units of the critical temper- 
ature T c ; the time is in units of the characteristic time 
to = 7rh/8KBT c [12[; k is the Ginzburg-Landau param- 
eter; r] and (3 are the relaxation times of the order pa- 
rameter and the vector potential, respectively. We have 
solved the TDGL equations upon using the link vari- 
ables method. 

[3 EH Since we consider invariance 
of the system along the z direction, our approach could 
only be applied to a square mesoscopic superconducting 
cylinder. However, it might also be used for a very thin 
superconducting film of thickness d < £(0) (d < 1 in di- 
mensionless units) , as long as we replace n 2 in Eq. [T] by 
K e ff = ft 2 /cL[l6[ The geometry we have considered is de- 
picted at the top left corner of Fig. [TJ the lateral size of 
the sample is ds, and of the antidot dAD- 

For the study presented in this Letter, the relaxation 
times were kept fixed at rj = f3 = 1. The GL parameter 
was also maintained at n e ff = 25, corresponding, for 
instance, to a very thin film of Nb with thickness d « 
7.3 nm (assuming £(0) = 40 nm, k = 2.125, n e ff = 

K 2 /(d/m) « 25). 

For T = 0.53, Fig. [T] shows a short hysteresis loop 
which was made for a mesoscopic superconducting square 
(hereafter referred to as sample Si) of dimensions ds = 
12, with a small antidot at the center, with size dAD = 2. 
Points marked with letters, (a), . . . , (/) , indicate the val- 
ues of H where a vortex (antivortex) either enters or exits 
the sample. The magnetic field is increased from zero un- 
til a value somewhat above that at which the first two 
vortices penetrate in (a). It is then reversed until the 
opposite value is achieved, and then reversed once more. 
At point (d), two antivortices penetrate the supercon- 
ductor. Points (b) and (e) correspond to the exit of a 
vortex and an antivortex, respectively, immediately after 
which one flux quantum still remains in the antidot. The 
point we will be mainly focusing on is (c) [or, equiva- 
lently, (/)]. This point corresponds to the entrance of an 
antivortex, which will encounter a vortex trapped in the 
antidot. The V-AV pair will annihilate either inside or 
outside the antidot. 

We have made a systematic study of the V-AV an- 
nihilation process, varying the experimentally accessible 
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FIG. 1: The magnetization curve as a function of the external 
applied magnetic field normalized to its maximum value. The 
meaning of the points marked in the figure are explained in 
the text. The arrows indicate the direction of the hysteresis 
loop. The inset illustrates the geometry of the problem. 

relevant parameters: T, ds, and dAD- The other param- 
eters were kept fixed as specified above. 

We turn now into a brief discussion of our results. 
Firstly, we found that for ds and dAD fixed, the V- 
AV annihilation occurs at a limited temperature inter- 
val. Above the upper limit of such interval, the vortex 
trapped at the antidot leaves the sample before the en- 
trance of the antivortex, and no collision can possibly oc- 
cur. For the sample Si this upper threshold is somewhat 
above T = 0.53. As a matter of fact, since the pene- 
tration length increases as the temperature approaches 
T c , the supercurrent density associated with the pinned 
vortex spreads to distances large enough to reach the ex- 
ternal border of the sample. This enhances an attractive 
interaction between the surface and the trapped vortex, 
facilitating its escape. 

At a given temperature, there is a minimum width 
of the superconducting frame to allow for a V-AV colli- 
sion outside the antidot. For example, for T = 0.23 and 
dAD = 2 the minimum size of the frame [13] is 6 (these 
dimensions define sample S2: ds = 8 and dAD — 2). Our 
reasoning is as follows: the antivortex nucleated at the 
external border is kept pinned by the surface barrier at 
the external border while, in turn, the vortex remains 
pinned at the antidot. Both entities attract each other 
and the interaction is larger for smaller superconducting 
frames. It is then easier to pull the vortex out of the 
antidot for shorter distances. On the other hand, if the 
superconductor width is increased, the antivortex starts 
its excursion toward the center before the V-AV attrac- 
tion becomes appreciable. Being thus less anchored by 
the surface, the antivortex travels all the way to the an- 
tidot and the annihilation process completes there. In 
support to this argument, Fig. [2] depicts the values of 
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the x component of the supercurrent density at the sur- 
faces y = 12 and y = 8 for both samples Si and £2, 
respectively. These values of J sx were calculated in the 
absence of applied field, H = 0, and thus represent the 
supercurrent due only to the vortex inside the antidot. 
The V-AV attraction (Lorentz force) will be proportional 
to J sx . Notice that, for both samples, the absolute value 
of J sx is maximum at the center of the square edge. One 
can clearly see that, in the middle of the sample edge, 
x = ds/2, the intensity of the Lorentz force is systemat- 
ically larger for the sample S2 than for Si. 

Fig. [3] shows the dynamics of the V-AV collision for 
the sample Si at T = 0.53. Immediately after entering 
the sample, only the antivortex is in the superconducting 
sea. As it moves forward to the antidot, at some instant 
the trapped vortex will come out and both will collide 
and annihilate. 
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FIG. 2: The supercurrent density along the top surface of 
the superconductor. The parameters used are specified in the 
figure. The inset shows the Lorentz force on the antivortex 
due to the vortex in the antidot. 

Another interesting characteristic we have found was 
that as the antivortex penetrates, its interaction with 
the vortex is so strong that both appear entirely de- 
formed. Usually, vortices look like a round flux filament 
surrounded by whirling super currents. However, when 
the vortex and the antivortex approach one another while 
still partially pinned at the corresponding surfaces, they 
become alongated, as if squashed, and form a narrow 
channel between the external surface and the antidot. 
This can be seen in Fig. [H which represents the topology 
of the order parameter i/j and the supercurrent distribu- 
tion at the very moment when the collision takes place 
for the sample Si at T = 0.53. It is very suggestive 
that, during the collision time, super currents nearly van- 
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FIG. 3: (Color online) The pictures show the topological 
shape of the phase of the order parameter; the phase goes 
from —TV (blue) to ty (red). The last picture (IV) corresponds 
to the stationary state of point (c) of the hysteresis loop of 
Fig.[ll 



ish accross the line of the elongated pair V-AV, just like 
in a phase-slip event [18], which is the dual process of 
Josephson tunneling. 




FIG. 4: (Color online) The absolute vale of the order param- 
eter (left) and the supercurrent distribution (right) for the 
sample Si; very dark color means order parameter very de- 
preciated. The size of the arrows (right) are not real; they 
were enlarged for better viewing. 

We also have estimated the instantaneous velocity of 
the vortex on its way from nucleation at the surface 
through its capture by the antidot (point (a) of Fig. [1]). 
In addition, the velocities of both the antivortex entering 
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FIG. 5: The modulus of the velocity of the vortex and the 
antivortex for the sample Si at T = 0.53 as a function of 
time (point (c) of Fig. [1}; the inset corresponds to point (a). 



average velocity slightly larger then that for a normal 
Abrikosov vortex: At = 17.2792t , Ay = 4.125f(0), 
which yields vy e = 6.8 x 10 3 m/s. A value somewhat 
larger - although in the same order of magnitude - might 
be attributed to the additional attraction exerted by the 
antidot. 

In summary, we have studied the dynamics of the a 
V-AV annihilation. The process was activated magneti- 
cally through a short hysteresis magnetization loop. Our 
results indicate that the V-AV collision can take place 
either inside the antidot or at the superconducting sea. 
In addition, we have shown that, during the annihilation 
process, vortices and antivortices may be stretched by 
a competition between their mutual attraction and the 
forces that act to keep them anchored at the surfaces. To 
complete the picture, we have also reported velocities of 
the entities while travelling through the superconducting 
frame of the sample. 

The authors thank the Brazilian Agencies FAPESP, 
FACEPE, and CNPq for financial support. 



the sample and that of the vortex exiting the antidot, on 
their route to annihilation (point (c) of Fig. [I]), were also 
determined. The results are shown in Fig.O It is clearly 
seen that the velocities achieved during the collision pro- 
cess are much larger than those of a vortex penetrating 
the sample (at least an order of magnitude, as shown in 
the inset). 

In Reference [l9| (see also [lo|, it was estimated 
that a phase-slip line can achieve a velocity which is 
two orders of magnitude higher than an Abrikosov vor- 
tex which, in turn, is one order of magnitude faster than 
Josephson vortices. As Abrikosov vortices are quoted [l9[ 
to travel at 10 3 m/s, this means that typical velocities for 
a phase-slip line are of the order of 10 5 m/s. Although 
there is a conceptual difference between the narrow chan- 
nel of Fig. [H (a composition of a vortex and an antivortex) 
and a phase-slip line, the similarity in format is reinforced 
by the fact that their velocities are approximately the 
same. For instance, for the sample 5i, the window dur- 
ing which the antivortex remains visible is At = 0.5517to 
and the distance it travels is Ay = 2.875^(0); inserting 
T c = 3.72 K and £(0) = 230 nm (the relevant param- 
eters for Sn, which were used by Sivakov et at. [191]), 
the average velocity is vav = 1.5 x 10 5 m/s. On the 
other hand, for the vortex exiting the antidot we ob- 
tain At = 0.2135£o and Ay = 1.625f(0); which gives 
vy s = 2.2 x 10 5 m/s. Thus, in terms of format and veloc- 
ity, the annihilation process could be envisaged as a phase 
slip. It is worth noticing that the large velocities antici- 
pated for the V-AV pair along the collision process, are 
similar to those developed during the early stage of a vor- 
tex avalanche, as observed by the authors of Ref. [6], who 
reported velocities as large as 1.8 x 10 5 m/s for avalanches 
in YBCO films. 

For the vortex first entering the antidot we found an 



[1] C. A. Duran, P. L. Gammel, R. E. Miller, and D. J. 
Bishop, Phys. Rev. B 52, 75 (1995). 

[2] P. Leiderer, J. Boneberg, P. Brull, V. Bujok, and S. Her- 
minghaus, Phys. Rev. Lett. 71, 2646 (1993). 

[3] T. H. Johansen, M. Baziljevich, D. V. Shantsev, P. E. 
Goa, Y. M. Galperin, H. N. Kang, H. J. Kim, E. M. Choi, 
M. S. Kim, et a/., Superconductor Science & Technology 
14, 726 (2001). 

[4] F. Colauto, E. M. Choi, J. Y. Lee, S. I. Lee, V. V. 
Yurchenko, T. H. Johansen, and W. A. Ortiz, Super- 
conductor Science & Technology 20, L48 (2007). 

[5] F. Colauto, E. J. Patino, M. G. Blamire, and W. O. 
Ortiz, Superconductor Science & Technology 21, 045018 
(2008). 

[6] B. Biehler, B.U. Runge, P. Leiderer, and R. G. Mints, 
Phys. Rev. B 72, 024532 (2005); B.U. Runge, U. Bolz, J. 
Eisenmenger, and P. Leiderer, Physica C 341-348, 2029 
(2000). 

[7] G. R. Berdiyorov, M. V. Milosevic, and F. M. Peeters, 

Phys. Rev. Lett. 96, 207001 (2006). 
[8] R. Geurts, M. V. Milosevic, and F. M. Peeters, Phys. 

Rev. Lett. 97, 137002 (2006). 
[9] V. R. Misko, V. M. Fomin, J. T. Devreese, and V. V. 

Moshchalkov, Phys. Rev. Lett. 90, 147003 (2003). 
[10] A. S. Melnikov, I. M. Nefedov, D. A. Ryzhov, I. A. Shere- 

shevskii, V. M. Vinokur, and P. P. Vysheslavtsev, Phys. 

Rev. B 65, 140503(R) (2002). 
[11] L. F. Chibotaru, A. Ceulemans, V. Bruyndoncx, and V. 

V. Moshchalkov, Nature (London) 408, 833 (2000). 
[12] N. R. Weethamer in Superconductivity, edited by R. D. 

Parks (New York, 1969). 
[13] W. D. Groop. H. G. Kaper, G. K. Leaf, D. M. Levine, 

M. Palumbo, and V. M. Vinokur, J. Comput. Phys. 123, 

254 (1996). 

[14] G. C. Buscaglia, C. Bolech, and A. Lopez in Connectivity 
and Superconductivity, edited by J. Berger e J. Rubin- 
stein (Springer, Heidelberg, 2000). 



5 



[15] E. Sardella, A. L. Malvezzi, P. N. Lisboa-Filho, and W. 

A. Ortiz, Phys. Rev. B 74, 014512 (2006). 
[16] D. J. Priour and H. A. Fertig, Phys. Rev. B 67, 054504 

(2003). 

[17] While changing sample and antidot sizes, we have always 

adopted unitary steps. 
[18] J. E. Mooij, and Yu. V. Nazarov, Nature Physics 2, 169 

(2006). 



[19] A. G. Sivakov, A. M. Glukhov, A. N. Omelyanchouk, Y. 

Koval, P. Miiller, and A. V. Ustinov, Phys. Rev. Lett. 

91, 267001 (2003). 
[20] A. Andronov, I. Gordion, V. Kurin, I. Nefedov, and I. 

Shereshevsky, Physica C 213C, 193 (1993). 
[21] S. G. Doettinger, S. Kittelberger, R. P. Huebener, and 

C. C. Tsuei, Phys. Rev. B 56, 14157 (1997). 



